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Abstract 

Wc propose a primordial-origin model for the composite configurations of global magnetic fields in spiral 
galaxies. We show that uniform tilted magnetic field wound up into a rotating disk galaxy can evolve into 
composite magnetic configurations comprising bisymmetric spiral (S=BSS), axisymmetric spiral (A=ASS), 
plane-reversed spiral (PR), and/or ring (R) fields in the disk, and vertical (V) fields in the center. By MHD 
simulations we show that these composite galactic fields are indeed created from weak primordial uniform 
field, and that the different configurations can co-exist in the same galaxy. We show that spiral fields 
trigger the growth of two-armed gaseous arms. The centrally accumulated vertical fields are twisted and 
produce jet toward the halo. We find that the more vertical was the initial uniform field, the stronger is 
the formed magnetic field in the galactic disk. 

Key v^rords: interstellar matter — magnetic fields; simulations — MHD; galaxies — spiral; galaxies 
— the Galaxy 



1. Introduction 

Global magnetic fields in disk galaxies are observed to 
show either bisymmetric spiral (ESS, hereafter S), ring 
(R) or axisymmetric (ASS, hereafter A) configuration 
(e.g., Sofue et al 1987). Figure 1 illustrates the observed 
magnetic configurations. Large-scale S fields on the order 
of '--^ 1 are observed in many spiral galaxies includ- 
ing the Milky Way (Sofue et al. 1986; Sofue et al. 1983; 
Han et al. 1998; Krause et al. 1989, 1998). An R field 
is observed in M31 (Berkhuijsen et al. 2003), while some 
galaxies such as IC 342 show A fields (Krause et al. 1989). 
In addition to the S configuration, the local field in the 
Milky Way are observed to show reversal of field direction 
with respect to the galactic plane from positive to nega- 
tive galactic latitudes (Sofue and Fujimoto 1983; Taylor et 
al. 2009). We call such configuration the galactic-plane- 
reversed (GPR) field. 

Edge-on galaxies often indicate V fields in their halos 
(Brandenburg et al. 1993; Krause et al. 2008). The 
central regions of spiral galaxies exhibit V fields as ob- 
served for M31 (Berkhuijsen et al. 2003) and M81 (Krause 
et al. 1989). The Milky Way's center indicates promi- 
nent V fields as strong as ^ 1 mG (Yusef-Zadeh et al. 
1984; Tsuboi et al. 1986; Sofue et al. 1986; LaRosa et 
al. 2005), where the magnetic energy density (pressure) 
much exceeds the thermal and turbulent energy densities 
of the interstellar gas, showing that the energy-density 
equipartition does not hold. The field strength is main- 
tained steady due to amplification by dynamo mechanism 
(Parker 1979) and annihilation by magnetic diffusion (e.g. 
Fujimoto and Sawa 1987). 



We emphasize that the large-scale magnetic configura- 
tions in spiral galaxies are a mixture of different types 
of configurations, and that they co-exist in one galaxy. 
Figure 1 illustrates the global field configuration in the 
Galaxy, which may generally apply to spiral galaxies. The 
origin of such composite topology has remained as a mys- 
tery. In this paper we try to solve this problem, and to 
explain the observed configurations as the fossil of primor- 
dial magnetic field when the galaxies were formed. We 
unify the basic topologies - S, A, RA, R and V - into a 
common origin. 

Magnetohydrodynamic (MHD) simulations of accretion 
disks have been extensively performed to study the dy- 
namics and jet formation in magnetized galactic disks 
and AGN torus around massive black holes (Kudoh et al. 
1998; Matsumoto 1999; Machida et al. 2000; Nishikori et 
al. 2006). Nishikori et al. (2006) traced the three dimen- 
sional MHD evolution of a galactic disk with a weak az- 
imuthal field, showing that the fields are amplified by the 
differential rotation as well as by the magneto-rotational 
instabilities (MRI: Balbus and Hawley 1991; Hawley and 
Balbus 1991). They showed that a tightly wound spi- 
ral configuration superposed by random fields is formed. 
Wang and Abel (2009) have recently performed MHD sim- 
ulations of formation of a disk galaxy and magnetization 
of interstellar medium. They have shown that a weak 
uniform field perpendicular to the disk is amplified to cre- 
ate micro G interstellar fields during the disk formation, 
which affects the structure formation of ISM and star for- 
mation. 

In this paper, we also try to confirm that our pro- 
posed scenario for the global magnetic configurations in- 
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Fig. 1. Schematic illustration of magnetic fields in spiral 
galaxies, including our Galaxy. 

deed works in a disk galaxy by three dimensional MHD 
simulations. We focus on the global topology of magnetic 
lines of force, and show that a unified initial field can 
evolve into various magnetic configurations as observed in 
spiral galaxies. The present paper would be, therefore, 
complimentary to the above cited papers, in which more 
detailed ISM physics and jet formations are discussed. 

2. Primordial Origin Hypothesis 

2.1. S (BSS) Field 

It has been suggested that the S configuration may be 
created by a rotating disk galaxy, if a larger scale primor- 
dial magnetic field of cosmological origin were wound up 
into the rotating disk (Sofue et al. 1986, 1987). Figure 2 
illustrates this idea: A large scale uniform magnetic field 
parallel to the galactic plane is wound up with the rotating 
disk. When the disk formation is completed, the S field 
is maintained with a constant pitch angle by magnetic re- 
connection near the neutral sheets. The field reconnection 
occurs in a time scale t determined hy t = A/Va, which we 
call the Alfven time. Here the reversal scale A of the field 
may be taken as the turbulent cloud scale between neigh- 
boring spiral arms on the order of ~ 100 pc. The Alfven 
speed Va is on the order of ~ 6 km s~^for B ^ 3/iG and 
\ atoms cm^"^. Then, we have i ^ 2 x 10^ years, which 
is sufficiently shorter than the rotation period of the ga 
laxy. 

2.2. R (Ring) Field 

It may happen that the primordial magnetic field before 
galaxy formation was not uniform or lopsided such as due 
to proper motion of the primordial galaxy. Then, a part 
of the amplified S field will be reconnected in the galactic 
disk, and creates an R field inside a certain radius, as 
illustrated in Fig. 2. Once a ring field is created, it will 
shrink to attain smaller radius due to the magnetic tension 
and disk gas accretion toward the center. Accordingly, the 
bisymmetric component, left over in the outer disk, will be 




Fig. 2. Cosmological origin of S field from a uniform field, 
and of R field and its amplification. 
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Fig. 3. Cosmological origin of V, A, and PR fields. 

wound again by the differential rotation of tlie galaxy, and 
another reconnection takes place, creating a second ring. 
Thus, the field strength of the ring component increases, 
as the galaxy rotates. 

2.3. A (ASS) Field 

Some galaxies show axi-symmetric (A) configuration, 
where the field direction is always away from or toward the 
center, as if there exists a magnetic monopole in the cen- 
ter. Such field configuration is interpreted as due to locally 
enhanced polarized emission inside the disk by sheared 
magnetic loops, each of which is closed between the disk 
and halo. A more simple explanation of the A configu- 
ration is a primordial-origin from a vertical uniform field 
as shown in Fig. 3. If a vertical primordial field was 
wound up by the galactic disk, it is twisted tightly near 
the galactic disk. The upper parts of the wound field with 
respect to the galactic plane have spiral configuration with 
their directions away from or toward the galactic center, 
and in the lower half of the disk the field directions are 
reversed. Since the radio emission is strongest near the 



galactic plane, the Faraday rotation takes place for the 
emissions transferring the upper half of the magnetized 
disk, and hence, the observed rotation measure (RM) ap- 
parently indicates A configuration. 

2.4- PR (Galactic Plane-Reversed) Field 

The distribution on the sky of Faraday rotation measure 
(RM) of quasars is often used to diagonyse the magnetic 
field configuration in the Milky Way (Sofue and Fujimoto 
1983; Taylor et al. 2009). The observed longitudinal vari- 
ation of RM, reversing from one arm direction to the next 
arm, is well explained by S type field. The RM distribu- 
tion also shows large-scale latitudinal variation, showing 
large-scale reversal of field direction from positive to neg- 
ative longitudes in the local arm, e.g. in the direction 
of I - 90° (Sofue and Fujimot 1983; Taylor et al. 2009). 
This indicates that the field direction is reversing in the 
upper and lower regions of the disk plane. Such galactic 
plane-reversed (PR) field may be related to the A field, as 
discussed below. The PR configuration may be created at 
the same time when an A field is created from a V field, 
as illustrated in figure 3. 

2.5. V (Vertical) Field 

The strong vertical field in the Galactic Center has been 
a mystery. One possible explanation is to attribute it to 
the fossil of primordial vertical field trapped to the galac- 
tic disk (Fig. 3). Suppose that a forming galactic disk 
was penetrated by a uniform magnetic field perpendicu- 
lar to the disk plane. According to the contraction of the 
galactic gas disk, the field is accumulated to the central 
disk. The field lines cannot escape across the disk: The 
diffusion time due to anomalous diffusion defined by the 
Alfvcn time calculated for a length scale A ^ i? ~ 3 kpc 
for a disk galaxy is ^ 5 x 10^ years, much larger than the 
rotation period. The vertical field is twisted by the disk 
rotation, and causes angular momentum loss, which accel- 
erates the disk contraction and amplify the V field. Thus, 
a strong V field is created in the nuclear gas disk, where 
the magnetic energy density becomes comparable to the 
kinetic energy density of gas due to the g alactic rotation. 
In fact the V field of a few mG in the Galactic Center has 
energy density as high as ^ 10~^ erg cm~'^, which is 
close to the kinetic energy density of the molecular gas by 
galactic rotation ^ 100 pc of Ey^ ~ 10^^ erg cm~'^, and is 
much larger than the thermal energy density Et ^ 10^^° 
erg cm^'^. 

2.6. Unified Scenario for Primordial Origin 

We, here, encounter a serious question why some or all 
of the different topologies, S (BSS), A (ASS), R (ring), PR 
(plane-reversed), and V (vertical), are observed simulta- 
neously in one galaxy. In order to answer this question, 
we propose a unified scenario of primordial magnetic field 
origin during the disk formation of a galaxy. We con- 
sider a tilted original magnetic field penetrating a form- 
ing galactic disk, as illustrated in figure 4. As the lines 
of force are wound up by the contracting gas disk, the 
strength of vertical component is amplified simply obey- 



4 



ing By ^ Boy{R/Ro)^'^ , where R is the radius and Rq is 
the original radius within which the original field was en- 
closed. When By reaches a critical strength, the tension 
of the field will become strong enough so that the toroidal 
(spiral) field is slipped off from the disk plane to catch up 
the original tilted off-plane field. The amplified vertical 
field is further accumulat ed to produce a strong V field 
in the central disk. 

Thus, the outer field is wound up to form an S field in 
the outer disk, and a V field in the central region. There 
appears a composite field orientation: V in the center, 
and S field in the disk. If there was some anisotropy su- 
perposed on the tilted uniform field, a part of the S com- 
ponent in the disk will be reconnected to produce an R 
field. This will result in a composite magnetic structure 
composed of the three configurations: V in the center, R 
in the inner disk, and S field in the outer disk. The ob- 
served asymmetry of the rotation measures with respect 
to the galactic plane (Sofuc and Fujimoto 1983) may be 
explained by an A component superposed on these con- 
figurations. 

If the original field was perpendicular to the disk, or the 
tilt angle was large enough, the wound-up field evolves in 
to spiral configuration with its direction always inward or 
outward in one side of the galactic disk. It shows the same 
structure in the other side of the disk, but the field direc- 
tion is opposite, e.g. outward or inward. Thus, the field 
configuration is ASS type in one side, and the same but 
reversed in the other side of the disk. If the field config- 
uration in the solar vicinity has such structure, we may 
explain the large-scale field reversal observed in the RM 
distribution on the sky. Also, linear polarization obser- 
vations of radio emission from such a disk with a plane- 
reversed ASS field will be observed to show a global A 
field, since the Faraday rotation measure is more weighted 
by the ISM and field in the nearer side to the observer. 

3. MHD Simulation 

3.1. The Method 

In order to clarify if the here proposed primordial- 
origin scenario for the galactic magnetism indeed works 
during galaxy formation, we performed three-dimensional 
magneto-hydro dynamical (MHD) numerical simulations. 
The basic equations to be solved are as follows. 



| + V(PV) = 0, 



P 



dv 



jxB 



— = V X V X B 770J, 

ot c 



(1) 

(2) 

(3) 
(4) 



where p, P, cf), v, B, j = cy xB/47r, T, t/q and S 
are the density, pressure, gravitational potential, veloc- 
ity, magnetic field, current density, temperature, mag- 
netic diffusivity and specific entropy, respectively. We 
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Fig. 4. Origin of S, R and V fields from uniform (top) and 
lopsided fields (bottom). 

solve these equations using the modified Lax-Wcndroff nu- 
merical method with the cylindrical coordinate (r, 0, z) 
(Matsumoto 1999; Nishikori et al. 2006). We here re- 
port only the essential points of the simulations in order 
to confirm the reality of the present model. Detailed de- 
scription of the numerical simulations will be given in a 
separate paper. 

We consider a gas disk rotating in a fixed gravita- 
tional potential of Miyamoto and Nagai (1975), which ap- 
proximately represent the observed rotation curve of the 
Galaxy and is expressed by 



(r,z) 



[r2 + {a, + (z2 + 62)i/2}2]i/2 



(5) 



where the parameters and bi are given in table 1. 

As the initial condition wc assume that the disk is pen- 
etrated by a uniform intergalactic magnetic field. The in- 
clination of the field lines to the rotation axis (z direction), 
9, is taken to be 0° (perpendicular to disk), 45° , 80° and 
90° (parallel). The interstellar gas is put freely, rotating 
in the potential with the centrifugal force being balanced 
with the gravity. The gas density is assumed to be expo- 
nential in the radial direction with a scale radius 7.26 kpc. 



No. 



Table 1. Parameters for the Miyamoto-Nagai potential. 



log(density) 



Component 


i 


tti (kpc) 


bi (kpc) 


M(M®) 


Bulge 


1 


0.0 


0.495 


2.05 X 10^° 


Disk 


2 


7.258 


0.520 


2.547 X 10" 




Table 2. Units adopted in this paper. 



Physical quantity Symbol 



Numerical unit 



Length 
Velocity 
Time 
Density 
Magnetic field 



vq 
to 
Pa 
Bo 



1 kpc 

207 km s-i 



4.7 X 10^ 



1.6 X 10-24 ^ -.^-3 



g cm 



VPov^ 26 ^lG 



The inner and the outer edges of the disk are 0.4 kpc and 
10 kpc, respectively. The disk gas is given an equivalent 
temperature of 10'^ K and a density of ~ 1 H cm^'^, so 
that the scale thickness of the disk is _ff ~ 100 pc. The 
intcrgalactic gas is assumed to have temperature of 10^ 
K. Figure 5 shows the initial density profile in the r and 
z direction. 

The grid numbers in the directions of radius, azimuthal 
angle, and height from the galactic plane, were 301, 66 
and 337, respectively. The grid spacing's in the radius 
was 0.01 kpc up to 1 kpc, and it is stretched for r > 1 
kpc. The spacing's in height direction is 0.002 kpc till 
0.15 kpc, and it is stretched for z > 0.15 kpc. The maxi- 
mum grid spacing for r and z is O.lkpc. In the azimuthal 
direction, it is uniform of 0.098 radian. The magnetic 
Reynolds number was taken to be i?niO = voro/vo = 2000, 
where ijo is the magnetic diffusivity. The corresponding 
magnetic diffusivity is 770 — 3 x lO^^cm^s"^. where we as- 
sume the turbulent diffusivity (e.g., Parker 1979). The 
initial plasma /3 value at r = = 1 kpc in the galactic 
plane was 100. 

The initial uniform field strength was given to be ^ 0.6 
yU G. The unit time of this simulation is the inverse of 
angular velocity, to = tq/vo = 4.7 My, at tq = 1 kpc and 
vo = {GMo/roy/^ = 207 km s^ifor Mq = 10^°M^. In or- 
der to mimic the galactic disk formation, the accumulated 
central gas is assumed to be transformed to stars, so that 
the gas density within central 200pc region is kept always 
low enough. 

3.2. Results of MHD simulations 

The result of MHD simulations is shown in Fig. 6. As 
the gas disk rotates differentially, the magnetic lines of 
force are twisted and wound up with the disk. When the 



Fig. 5. Initial gas density profile in the {r,z) plane. 

field lines are wound tightly enough, they are reconnected 
at the neutral sheet within the magnetic diffusion time, 
which is on the order of ten times the crossing time of 
Alfvcn wave across the disk thickness d. The diffusion 
time is longer in the outer disk, while it is shorter in the 
central region, because the disk thickness is larger in the 
outer disk. 

Figures 8, 9 and 10 show the distributions of logarithm 
of the gas density and magnetic energy density for the case 
of = 0°, 45°, and 90°, respectively. The gas and magnetic 
energy densities are averaged around the equatorial plane 
(|z| <0.02). 

3.3. Evolution of vertical field 

In the case of 6* = 0°, the vertical field lines that pen- 
etrated the outer disk are wound up to produce spiral 
field whose direction reverses at the galactic plane. In the 
upper side of the disk, the field lines point inwards (or out- 
wards) at any azimuthal angle, showing A configuration, 
whose lower half shows opposite field direction. In the 
central region of the disk, the initial field lines are twisted 
to create a helical V field. The global field configuration 
becomes, thus, a composite of outer A and inner V fields 
as shown in Fig. 7. According to the twisting, the upper 
and lower inversely directing components are reconnected 
near the disk plane, which results in accumulation of the 
vertical field lines toward the central region. The accu- 
mulated vertical fields are further twisted in the center, 
creating stronger vertical field. 

Thus a very strong vertical field is created in the cen- 
tral region. In the present computations, the elapsed 
time corresponds to a few hundred million years, or one 
galactic rotation time near the solar orbit. However, the 
galactic rotation within the central 1 kpc is on the order 
of ten million years, so that the disk rotates almost ten 
times. In such rotation periods, the vertical component 
is strongly accumulated, and strength is amplified by a 
factor of (r/ro)2, where r is the final radius of the vertical 
field tube, and rg is its initial value. If tq ^ 10 kpc, and 
the r ~ 1 kpc, the initial field is amplified by a factor of 
10^. Hence, even if the initial field was very weak, very 
strong central vertical field is created. 

According to the accumulation and twisting, the field 
strength near the center is significantly amplified, and the 
magnetic pressure gradient perpendicular to the disk be- 
comes strong enough to accelerate the gas toward the halo. 
This causes vertical outfiow of gas as well as the loss of 
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Fig. 6. MHD simulation of the primordial origin model of magnetic field in a disk galaxy. Curves show the magnetic lines of force. 
Blue and green surface show the iso-surface of the density at p = 0.1 and 0.28, respectively, {a)-(c): Evolution of field lines originally 
incUned by 9 = 45° to the rotation axis at t = 0, 14.5 and 30.0 to with to = 4.7 My. (d)-(f): Field lines for 6» = 0°, 80° at 14.5to and 
90° at t = SO.Oto. Note that S or A, and V fields are formed in the same galaxy. 




Fig. 7. Central region of Fig. 6d. 
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1=19,0 e=0 log\/3) log(ByS7T) 




Fig. 8. Gas density log p (left) and magnetic energy density 
log B-^/Stt (right) for the case of = 0. Here, B and p are 
normalized by Bq = 26 fiG and po = 1 H cm"''. The gas 
density is maximum near the neutral sheet, where the field 
strength is minimum. 



1 = 30,0 = 45 log(^'} log(ByS7T) 




Fig. 9. Same as fig. 8, but for the case oi = 45°. 




Fig. 10. Same as fig. 8, but for the case oi 9 = 90°. 

angular momentum of the disk by magnetic braking, ac- 
celerating the contraction of the galactic disk (Uchida et 
al. 1986). 

3.4- Initially inclined field: V, A and PR fields 

The cases for 6 = 45° and 80° show that tilted lines of 
force are wound up to create S configuration near the disk 
plane. The field lines penetrating the inner disk are also 
wound up significantly and their tension force makes the 
magnetic field "stand up". The twisting magnetic field 
lines remove the angular momentum of the gas, and the 
gas is accreted to the central region with the field line. 
Eventually, due to the dissipation and upward propaga- 
tion of the tightly wound up toroidal magnetic field in the 
central region, V field is produced. 

In the case of 6* = 45°, the ratio of each component of the 
magnetic field strength is about \Br\ : |-B</,| : \Bz\ '-^ 3 : 10 : 1 
in the central region (0.3 < r < 0.8). The ratio is similar to 
that obtained in the accretion disk in which the magneto- 



rotational instability is developed. 

The highly inclined field with 6 = 80° similarly creates S 
in the disk and V field in the central region. Thus, outer 
S and inner V configurations of galactic magnetic fields 
are created simultaneously from a tilted primordial field, 
not strongly dependent on the inclination angle. 

When the initial field was more vertical, an A type con- 
figuration is dominant in the outer disk. It should be 
noted that the A field direction is reversed with respect 
to the galactic plane. Such plane-reversed (PR) field is 
indeed observed in the local disk of the Milky Way. 

3.5. Initially parallel field and Grand- designed spiral 
arms: V and S field 

The case for 9 = 90°, e.g. for a primordial field parallel 
to the disk plane, results in an S field almost over the en- 
tire disk. The field intensity map shows that two-armed 
grand-design spiral pattern is formed as shown in Fig 10, 
where B'^ /Sir is shown by grey scale intensity and is com- 
pared with the gas density. It must be stressed that the 
gaseous spiral arms are formed along the minimum loci 
of the field intensity (neutral sheet). The similar density 
structure is also appeared for various inclination angles 
except for 9 = 0° (see Fig. 8 and 9) In the central region, 
the field lines are tightly wound, resulting in high mag- 
netic pressure. This causes infiation of magnetic flux, and 
also creates the V field in the central region (Fig. 6(f)). 
The ratio of each component of the magnetic field strength 
is about |Br-| : l-B^I : 1 : 100 : 10 in the central region 

(0.3 < r < 0.8). The ratio shows that the field structure is 
almost R field in the central region. 

3. 6. Relation of created fields and the initial field incli- 
nation 

Figure 11 shows the time evolution of the magnetic en- 
ergy per volume, B'^ /Sir, for various inclination of the 
initial field direction. The magnetic energy density per 
volume is averaged within the central region of the disk 
(0.3 < r < 0.8, \z\ < 0.04). The field strength for initially 
vertical field is strongly amplified. On the other hand, if 
the initial field was parallel to the disk, the azimuthal field 
is not drastically amplified. 

The strength of the created spiral field is weaker than 
that for smaller inclination, because the parallel field lines 
to the disk plane can be easily dissipated when it is tightly 
wounded. The larger is the inclination (more vertical) of 
the original field, the stronger is the created spiral field. 
In other words, the centrally accumulated vertical field is 
proportional to the total flux of original fleld penetrating 
the initial galactic disk. If the disk thickness is ignored, 
the total flux F, and therefore, the final amplified field 
strength is roughly represented by ~ Focos9. Namely, 
the total flux was larger for larger inclination, and vice 
versa. 

This relation between the initial angle of the primor- 
dial field and the resulting field strength predicts that, if 
the initial field was perfectly parallel to the disk plane, 
e.g. 9 = 90°, the resultant field strength would be very 
weak. Also, the total flux of magnetic fleld of a galaxy 
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Fig. 11. Time evolution of the magnetic energy density per 
volume /Sir. Initially vertical fields are more strongly am- 
plified compared with initially parallel fields. 



Fig. 13. Accretion rate of the disk gas as functions of time. 
Accretion, and therefore the magnetic braking, is more effi- 
cient for initially vertical field compared to the initially par- 
allel field. . 
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Fig. 12. Time evolution of plasma /3 value, e.g. the ratio of 
the gaseous to magnetic energy densities. The more vertical 
is the initial field, the stronger is the decrease of /3. 

reflects the initial field inclination. We may thus intro- 
duce a field-strength function for spiral galaxies depend- 
ing on the initial inclination. If the initial disk direc- 
tions of primordial galaxies were randomly distributed 
with respect to the large-scale cosmic magnetic field, the 
number density of galaxies with total field flux F would 
be represented by the probability distribution of angle 
e as dN{F)/dF = Csin6l = C 1 ~ {F / Fo)^_. Instead of 
F we may use the mean field strength B and Bq as 
dN{B)/dB = C"y/l-B/Bo. 

3.7. Initial field direction and the magnetization of the 
ISM 

Figure 12 shows the time evolution of plasma (3, the 
ratio of the thermal energy density of gas to magnetic 
energy density in the disk. The plasma /3 is also averaged 
value within the central region of the disk (0.3 < r < 0.8, 
l^l < 0.04). The more vertical was the initial field to the 
disk, the more strongly decreases the (3. The reason is 
simply because of the initial magnetic flux penetrating 
the whole disk, as discussed in the previous subsection 
related to the azimuthal field amplification. 



3.8. Magnetic accretion of the disk, angular momentum 
transfer, and jet formation 

Figure 13 shows the time variation of accretion rate of 
gas in the disk for the various inclination of the initial 
field direction. The accretion of gas takes place due to 
the angular momentum transfer due to the magnetic ten- 
sion by twisted vertical field anchored to the intergalac- 
tic space through the j x B term in equation of motion. 
Accordingly, the accretion rate is strongest for the case of 
initial field perpendicular to the disk, as shown in figure 
13, while it is low for parallel field case. The accretion 
fiow will enhance the vertical component of the magnetic 
filed in the very central region of the galaxies (< 100 pc), 
although it is not resolved in the presenting simulations. 

The twisted field yields magnetic pressure gradient in 
the direction perpendicular to the disk, and produces a 
polar cosmic jet. The jet formation is strongest, when the 
initial field is vertical. This is a well known mechanism to 
create a cosmic jet since Uchida and Shibata (1985). It 
should be emphasized that even for a nearly parallel initial 
field with 9 ^ 80° the accreted vertical component in the 
galactic center is amplified enough to cause a polar jet. 
Thus, the vertical field and jet formation near the nuclei 
are common phenomenon in almost all spiral galaxies, if 
the present primordial origin hypothesis applies. 

3.9. Field amplification and magneto-rotational instabil- 
ity 

The magnetic Reynolds number in a rotating disk is de- 
fined by Sano & Miyama (1999) as i?,„c = v\^/{r]oV,), and 
is a good indicator to investigate the magneto-rotational 
instability (MRI). Here, vaz = ^zl ^^t^P-, ^ is angular 
frequency of the disk, ryo the magnetic diffusivity, and 
Bz the z component of the magnetic field. When 
is less than unity, the dissipation process generally sup- 
presses the growth of instability. The magnetic Reynolds 
number may be expressed in terms of i?niO as i?mc = 
i?mo(cs/wo)^(wA2/cs)^, whcrc Cs is the sound speed in the 
disk. Since the temperature in the disk is assumed to be 
~ lO'' K, we have Cs/wo ^ H/r^ ^ 0.05 in the disk, where H 
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Fig. 14. Time evolution of magnetic energy density for individual magnetic components. Thick solid, long-dashed and thin-dotted 
lines show _B^/(87r), B^/{8n), and B'^/{Sn), respectively. The four panels indicate the cases for different inclination angles: 9 = 0°, 
d = 45°, 9 = 80°, and 9 = 90°, respectively, from top left to bottom right. 



is the scale height of the disk. In the case of = 0° , we have 
(vAz/cs)'^ = {'^1 Pi) — 1/50 at the initial condition, which 
leads to i?mc — 0.1. The wavelength for maximum insta- 
bility is given by Am ~ ^/i'Az ^ O.OTrg. This wavelength is 
slightly larger than the scale height of the disk. Therefore, 
the axi-symmetric mode of the MRI is suppressed by the 
magnetic diffusivity in the present case. 

Even when the MRI is suppressed inside the disk, it 
can develop near the surface of the disk (Sano & Miyama 
1999). Near the disk surface, the wavelength for maximum 
instability becomes shorter than the scale height, because 
the Alfven speed is larger for the lower density. Once the 
MRI develops there, surface accretion occurs even in the 
resistive disk (Kuwabara et al. 2000, 2005). These pro- 
cesses amplify the strength of the magnetic field. Fig. 14 
shows the time evolution of the magnetic energy density 
for each component of the magnetic field, where the mag- 
netic energy densities are the values averaged within the 
disk at 0.3 < r < 0.8 and \z\ < 0.04. Each panel shows the 
result for different inclination angle of the initial magnetic 
field. 

In the case of 6* = 45°, the initial amplification of the 
magnetic field (around t < 5) is caused by the differen- 
tial rotation of the disk and the energy density of 
increases at first. It begins to increase almost exponen- 
tially around after t ^ 8 when the accretion has become 
prominent (see Fig. 13 ). The energy density of B,. in- 
creases almost exponentially from the beginning. The en- 
ergy density of also slightly increases, and the ratio of 



between the individual magnetic components remains at 
about \Br\ : \B^\ : \Bz\ ~' 3 : 10 : 1. The ratio is a typical 
value that is obtained in an accretion disk with develop- 
ing MRI (Balbus and Hawley 1991). We also found that 
the speed of the accretion fiow near the disk surface is 
larger than that at 2 = where the flow speed is some- 
times positive in the r direction. Thus, we may conclude 
that the magneto-rotational instability develops near the 
surface of the disk. The total magnetic energy is eventu- 
ally saturated around t ~ 20. In the case of 9 — 0° , the 
time evolution is similar to that of 6* = 45°, although the 
magnetic energies are not saturated yet. 

The growth rate of the field strength, as represented by 
the time in which the amplitude increases by a factor of e, 
is found to be approximately equal to 0.1/to from figure 14 
and figure 11. This rate may be compared with the growth 
rate of perturbation as predicted from the linear theory of 
MRI. The maximum local growth rate predicted from the 
rotation curve of the Miyamoto-Nagai potential is about 
0.7/io at 0.3 < r/ro < 1 within an ideal MHD approxima- 
tion. The growth rate as we found here is slightly smaller 
than this maximum value, probably because the MRI in 
our simulation is suppressed by the the resistive dissipa- 
tion process. The result is consistent with the global linear 
theory by Sano & Miyama (1999) which studied MRI in 
a resistive stratified Keplcrian disk with vertical magnetic 
field. The global linear analysis of the surface accretion 
has not been done yet for our initial setup. The more 
quantitative comparison with the global linear analysis 
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would be subject for future work. 

In the case of 6* = 90°, the initial amplification is also 
caused by the differential rotation of disk, as in the case 
of 6* = 45°. Then, the energy density of decreases 
around < ~ 10 because of the magnetic diffusion or by re- 
connection in the disk. The energy density of also de- 
creases because of the magnetic diffusion or reconnection. 
A small amount of B, field is amplified in the disk around 
t ^ 13 possibly because of the vertical motion driven by 
the thermal pressure during the magnetic diffusion or re- 
connection. The Bz field is stretched by the rotation of 
the disk, and B^ field is amplified again. Because the 
accretion flow does not develop in time (see Fig. 13), Br 
does not increase. In the case oi 9 = 90°, the MRI near 
the surface of the disk does not develop because there is 
no global vertical field crossing the disk. Each component 
of magnetic energy eventually saturates around i ^ 20 and 
keeps almost constant during our calculation. In the case 
of 6* = 80°, time evolution is something between that of 
6 = 45° and 6 = 90°, although the magnetic energies have 
not saturated yet. 

Even though the axi-symmetric mode of the MRI is 
suppressed in the disk, the non axi-symmetric mode of 
MRI may develop when a weak toroidal component of the 
magnetic field is dominated in the disk (e.g., Balbus and 
Hawley 1992; Terquem & Papaloizou 1996). However, 
there is no indication of non axi-symmetric mode of the 
MRI, even in the case oi9 = 90°. The most- unstable wave- 
length of the non axi-symmetric mode may not be resolved 
in our simulation, because only 64 grid points are used in 
the 4> direction. The longer wavelength may develop later, 
but it may take time because the linear growth rate of the 
non axi-symmetric MRI is generally smaller than that of 
the axi-symmetric one (e.g., Matsumoto & Tajima 1995). 
The further simulations with fine grid size and longer time 
scale will be needed in future. 

4. Discussion 

We have proposed a new scenario, which we call the 
primordial-origin model, in order to explain the observed 
composite magnetic configurations in spiral galaxies. A 
tilted weak uniform magnetic field is wound up by a ro- 
tating disk galaxy and creates the four basic magnetic 
topologies, S, A, R and V (bisymmetric spiral, axisym- 
metric spiral, ring, and vertical). We also showed that 
the reversal of field direction with respect to the galactic 
plane (PR - plane reversal - configuration), as observed 
in the local disk of the Galaxy, is formed from the verti- 
cal component of primordial field during the formation of 
A configuration. In order to show that this scenario in- 
deed applies, we have performed three-dimensional MHD 
simulations of magnetized rotating gas disk, mimicking a 
forming disk galaxy. 

We found that the more vertical was the initial field, 
the stronger is the created field. The V field is created 
and amplified in the central region even from tilted uni- 
form field. The magnetic fields infiuencc the interstellar 
gas, and the twisted vertical field accelerates outfiow of 



the central gas toward the halo. This leads to angular 
momentum loss of the disk, and enhances the disk accre- 
tion. The S field creates two-armed spiral arm structure, 
and the gas is compressed along the spiral neutral sheets, 
yielding a grand-designed spiral structure in the gaseous 
disk. 

We focused on the global topology of galactic magnetic 
fields, and explained the variety of observed configurations 
from a unified primordial-origin hypothesis. On the other 
hand, detailed analyses of the interstellar physics arc still 
crude. In this context our paper may be supplementary to 
the previous MHD simulations of magnetized disks such as 
those by Kudoh et al. (1998), Matsumoto (1999), Machida 
et al. (2000), Nishikori et al. (2006), and more recently 
by Wang and Abel (2009), in which detailed analyses and 
discussions are given for the ISM physics. 

The main purpose of this paper was to propose the idea 
for explaining the origin of the different global magnetic 
configurations in spiral galaxies. The idea was basically 
confirmed by MHD simulations of the global galactic disk 
with the empirical turbulent magnetic diffusivity. The 
numerical simulation showed that the vertical (V) struc- 
ture of magnetic field in the central region of the galaxy 
is created, even when the initial inclination angle of the 
magnetic field is nearly parallel to the galactic disk. It 
also showed that the S, A, or PR magnetic field topolo- 
gies are realized depending on the initial conditions. The 
magnetic energy density was found to be saturated after 
several rotations in the central region of the disk. The 
averaged amplitude of the saturated magnetic field in the 
central region (< Ikpc) was about ^ T^G in the case of 
6* = 45°, and - 2^G in the case of 6* = 90°. 

The amplification of the global magnetic field is caused 
by the differential rotation of the disk as well as the ac- 
cretion flow to the central region. The field strength gets 
saturated by the turbulent diffusivity after some rota- 
tions in the central disk. However, in order to discuss 
the strength of the magnetic field in the solar neighbor- 
hood, we need further simulations with longer time scales. 
The galactic dynamo process may be needed to get fur- 
ther amplification of the magnetic field (e.g., Hanasz et 
al. 2009). We have also considered possible growth of 
the magneto- rotatational instability (MRI). The instabil- 
ity is suppressed inside the galactic disk around the galac- 
tic plane. However, it was found that the MRI grows near 
the surface of the disk, and magnetic fields are also am- 
plified by MRI. 

In the MHD simulation, we have assumed an initially 
uniform magnetic field as the primordial seed field. This 
uniform assumption may apply in such cases that the uni- 
formity of the field was sufficiently larger than the scale 
radius of proto-galactic disk. Namely, the present sce- 
nario and the MHD simulation may be applied to such cir- 
cumstances that intergalactic magnetic fields existed with 
scale lengths comparable to the size of a cluster of galaxies. 
Since the major interest was to see the topological evolu- 
tion of initially uniform magnetic field in a single galaxy, 
the calculated area covered only a few tens of kpc around 
the galaxy. Moreover, the initial field was assumed to be 
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already as strong as ~ 0.6//G. In this context, the present 
simulation was not aimed at fully tracing the cosmological 
galaxy formation, starting from much smaller density per- 
turbations and weaker intergalactic field before structural 
formation in the universe. 

In order to see if a much weaker primordial magnetic 
field can be amplified to a micro G level during the galaxy 
formation, we here try an order-of-magnitude estimate of 
possible evolution of the field strength along the idea of 
primordial origin hypothesis (see also Sofue and Fujimoto 
(1987)). Suppose that an intergalactic primordial mag- 
netic field of one nG was trapped to a pre-galactic sphere 
of Tp ^ 100 kpc radius, and that the sphere collapsed 
to a proto-galactic disk of ~ 10 kpc. Then the field 
strength is amplified by a factor of (rd/rp)^ ~ 100, yielding 
0.1/xG. If a half of the disk flux was further accumulated 
to the central region, forming a proto bulge before star 
formation, of radius rb ~ 1 kpc, it is further amplified by 
(rd/rb)^ ^ 100 times. Thus, amplification by a factor of 
100 X 100 = 10"* times may easily be attained to create a 
central disk field of ^ lOfiG. Moreover, the twist by the 
galactic rotation, as well as accretion by magnetic brak- 
ing, would further amplify the vertical field by a factor of 
10 to 100 times, creating vertical fields of a few mG field. 
Thus, we may have ~ 10^ times amplification to a few fj,G 
in the disk, and further ~ 10® times amplification to mG 
field in the center, starting from intergalactic primordial 
field on the order of ~ InG. 

Finally we comment that the self gravity, star forma- 
tion, and more sophisticated ISM conditions such as the 
existence of multi-phase components and their phase tran- 
sitions in a more realistic galaxy-formation model should 
be subject for future works. The present simulation is, 
therefore, preliminary and should be taken as the initial- 
attack toward more general studies of the origin of galactic 
magnetism in the scheme of galaxy formation, considering 
the self-gravity and cosmological magnetic field origin. 

Numerical calculations were carried out on NEC SX-9 
at the Center for Computational Astrophysics (CfCA) of 
National Astronomical Observatory of Japan 
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